Solitons in two-dimensional Bose-Einstein condensates 
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The excitations of a two-dimensional (2D) Bose-Einstein condensate in the presence of a soliton 
are studied by solving the Kadomtsev-Petviashvili equation which is valid when the velocity of the 
soliton approaches the speed of sound. The excitation spectrum is found to contain states which 
are localized near the soliton and have a dispersion law similar to the one of the stable branch of 
transverse oscillations of a ID gray soliton in a 2D condensate. By using the stabilization method 
we show that these localized excitations behave as resonant states coupled to the continuum of free 
excitations of the condensate. 
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Bose-Einstein condensates of ultracold atoms are ideal 
systems for exploring matter wave solitons [1]. For 
most purposes, these condensates at zero temperature are 
well described by the Gross-Pitaevskii (GP) equation Q 
which has the form of a nonlinear Schrodinger equation, 
where nonlinearity comes from the interaction between 
atoms. In one dimension, the GP equation for conden- 
sates with repulsive interaction admits solitonic solutions 
corresponding to a local density depletion, namely gray 
and dark solitons. Such solitons have already been cre- 
ated and observed in elongated condensates with diverse 
techniques 

Also multidimensional solitons in condensates have at- 
tracted much attention [4j . An interesting type of excita- 
tion in a two-dimensional (2D) condensate is represented 
by a self-propelled vortex-antivortex pair which is apar- 
ticular solitonic solution of the GP equation 0, 0, 0, @] . 
In the low momentum limit, the relation between the en- 
ergy E and momentum P of the soliton [5[ approaches 
the dispersion law of Bogoliubov phonons, e = cq, from 
below. In this limit, when the 2D soliton moves at a ve- 
locity V close to the Bogoliubov sound speed c, the phase 
singularities of the vortex-antivortex pair disappear and 
the soliton takes the form of a localized density deple- 
tion, also called rarefaction pulse. Moreover, if V is close 
to c the GP equation can be rewritten in a simpler form, 
known as the Kadomtsev-Petviashvili (KP) equation 

In a previous short paper [10], we already presented 
some preliminary results on the dynamics and stability 
of this 2D soliton. In particular, we studied the excita- 
tions of the condensate in the presence of the soliton by 
linearizing the KP equation around the stationary solu- 
tion. By looking at the shape of the eigenfunctions we 
found excitations localized near the soliton, having shape 
and dispersion law similar to those of the transverse os- 
cillations of a ID gray soliton in a 2D condensate. In 
this work we present a more systematic analysis. We use 
a stabilization method in order to obtain a better de- 



termination of the dispersion law of the localized states. 
Moreover, the same method allows us to visualize the 
coupling between the localized states and the free states, 
i.e., the Bogoliubov phonons of the uniform condensate. 

Let us first summarize the derivation of the KP equa- 
tion. We consider a 2D condensate with a soliton moving 
at a constant velocity V in the x-direction. If the den- 
sity at large distances is noo, one can define the healing 
length £ = fi/[2m^n 00 ] 1 / 2 , where g is the mean- field cou- 
pling constant and m is the mass of the bosons. One can 
also introduce the dimensionless variables x — > y — > 
and t — > mt£ 2 /h, the normalized order parameter 
* -> V^o^> and tne velocity U m^V/h V/(cy/2). 
In the frame moving with the soliton, the GP equation 
isH 
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The order parameter can be written in the form \£ 



1/2 



exp[zS]. The equations for the phase and the density 
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When V is close to c, one has U 
a small parameter e = vT 

1/(2 V^)^ 2 , and expand the density and the phase in this 
form 



1 / y/2. Let us introduce 
2Z7 2 , so that U ~ l/y/2- 
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To the lowest order in e, the GP equation gives ds/dx — 
-ef/y/2 and 
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where we have introduced the stretched variables x = 
-ex + eH/(2y/2), y = e 2 y/V2 and t = e 3 t/(4 v / 2). Equa- 
tion ([6]) is known as the Kadomtsev-Petviashvili (KP) 
equation Hereafter we will omit the tilde in the 

stretched variables. 

The KP equation (j6j) admits a stationary solution of 
the form 



fo(x - 2t) = sech 2 \(x - 2t)/V2 



(7) 



which is independent of y and is a limiting form of a 
ID gray soliton The linear stability of this gray 

soliton in two dimensions was studied in 

[la El- One 

can look for transverse fluctuations propagating along y 
of the form f(x,y,t) = f (x - 2t) + ^{x - 2t)e { ^-^ 
and solve the KP equation up to terms linear in ip. One 
obtains the dispersion law 
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FIG. 1: Solitonic solution @ of the KP equation © as a 
function of the stretched dimensionless lengths x and y (see 
text for definitions) . For a given value of the small parameter 
e, the density distribution of the condensate is given by n — 
1 — e 2 f and the ratio between the soliton widths along x and 
y, in dimensional units, scales as e. 



The frequency uo is real for k > y/3/2 and imaginary for 
k < y/S/2. The gray soliton ([7j) is thus unstable: long 
wavelength transverse oscillations can grow exponentially 
(snake instability). We note also that for k ^> 1 the stable 
branch of excitations behaves as uj ex & 3 / 2 , similarly to 
the dispersion of capillary waves on the surface of a liquid. 

The KP equation (|6|) admits also a 2D solitonic solu- 
tion of the form lljl 



fo(x - 2t,y) 



4[f + 2y 2 -jx- 2tf 
[§ + 2y 2 + (x - 2t) 
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(9) 



which corresponds to the V — > c limit of the rarefaction 
pulse The function /q is plotted in Fig. 1. Differently 
from the ID gray soliton 0, the solution ([9j) decays al- 
gebraically in all directions. 

In the units of the GP equation ([T]), the density profile 
of the soliton is found by inserting Eq. ([9|) into Eq. (|4]). 
Its size along x and y is of the order of 1/e and l/e 2 , 
respectively. To the first order in e, the momentum P 
and the energy E per unit length are Q 
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As it must be, the soliton has a sound-like dispersion in 
the limit where the KP equation is valid, that is, when 

Now we consider density fluctuations of the form 
f(x,y,t) = f (x - 2t,y) + i/j(x - 2t,y)e~ iujt . Inserting 
this expression into Eq. ([6]) and keeping terms linear in 
ip, we find 
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It is convenient to introduce the function y) such 
that %p = d 2 (j)/dx 2 . By integrating Eq. JTTJ) twice in x, 
one has 
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We numerically solve this equation in a box of size 
L x x L y . The function <p is expanded in plane 
waves with periodic boundary conditions: <p(x, y) = 

^^xiA^xiAv), with X iAA = L~ 1/2 e i2 ™ x i L * 

for \v\ < l x (y 7^ 0). Concerning X2,/j,(y) we note 
that Eqs. (fTTj) and ([T2|) are invariant for y — > —y. 
Therefore the function <p is either an even or odd func- 
tion of y. If it is even, then one can take X2,/j,(y) — 
(2/Ly) 1 / 2 cos(27Tfiy/L y ) for 1 < fi < l y and X2,o(y) = 
(1/Ly) 1 / 2 . If it is odd, one can take X2,n(y) — 
(2jLy) x l 2 sin(27r jiiy / L) for 1 < \i < I. One thus obtains 
the following matrix equation 

i-ql - 2q x - ^ + 6 J2 —M^^^w* 

= LO(p v ^ (13) 

/Las/2 rLy/2 
dx / dy Xi,i/X2,/x/oXi,i/ , X2,^,(14) 
-L x /2 J-L y /2 



where q x = 2iiv/L Xl q y = 2nn/L. The size of the matrix 
is Nx TV, where N = 2l x (l y + l) for <p even, and N = 2l x l y 
for (j) odd. 

By solving Eq. (p~3]) . we find that all eigenvalues are 
real and positive. This is consistent with the stab ility 
of the 2D soliton 0, EEf. As already discussed in [Toj . 
among the eigenvectors we find states which are local- 
ized near the soliton. In order to better identify these 
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FIG. 2: Eigenfrequencies u> of Eq. (|13[) as a function of L y for 
L x = 10. The lower plot is a magnification of the upper one 
around u) = 4. All quantities are dimensionless. 

states and explore their coupling with the free, unbound 
states we use a stabilization method [l6|. In this method, 
the eigenvalues are calculated repeatedly for different box 
sizes L y and L x . Once the frequency is plotted as a func- 
tion of the box size, the localized states are identified as 
those having a dispersion which becomes fiat for large 
boxes, larger than the typical size of the corresponding 
bound state. These states are immersed into a bath of 
unbound states. The latter are characterized by the num- 
ber of nodes of the eigenfunctions in the x and y direc- 
tions and, due to the finite box, they appear as series 
of discrete branches in the stabilization diagram. The 
dispersion laws of the unbound states, as well as their 
dependence on L y and L x , can be estimated analytically 
by inserting fo = into the equations (p~3|) and A 
coupling between bound (resonant) states and unbound 
(free) states can be seen in the stabilization diagram in 
the form of avoided crossings due to non-zero matrix el- 
ements connecting the two types of states [l6| . A typical 
diagram is shown in Fig. [2j where we plot the eigenfre- 
quencies of Eq. ([T3]) as a function of L y for L x = 10. Flat 
dispersions and avoided crossing are clearly visible. An 
example of localized state is shown in Fig. 03 where we 
plot the x = section of the wave function of one of the 
bound states of Fig. [2] together with the soliton profile. 

In the soliton region the wave function of the excited 
states exhibits oscillations in the y direction whose wave 
vector k can be easily estimated. It is interesting to plot 
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FIG. 3: Wave function ip(0,y) (solid line) of an excited state 
localized near the soliton. The soliton profile, — /o(0,?/), is 
shown as a dashed line. The excited state has frequency oo = 
4.05 and the box size is L x = 10 and L y = 23. All quantities 
are dimensionless. 

the eigenfrequency uo of the bound states, extracted from 
the stabilization diagram, as a function of /c, extracted 
from the shape of the corresponding wave functions. This 
is done in Fig. [4] for the bound states having the lowest 
number of nodes in the x directions. Our results are 
shown as points with error bars, where the error bars are 
of the order of the inverse of the size of the soliton along 
y. In the same figure, the dispersion law (|HJ) of the stable 
branch of excitations of a ID gray soliton is plotted. The 
spectrum is very similar. This reflects the fact that a 2D 
soliton moving with velocity close to c is very elongated 
in the y direction (the ratio between the widths along x 
and y, in dimensional units, is proportional to the small 
parameter e) and its density distribution is indeed sim- 
ilar to that of a ID gray soliton. Differently from a ID 
gray soliton the 2D soliton has a discrete spectrum of 
bound states due to its finite length. This finite length 
also implies an "infrared cutoff in the spectrum of the 
bound states, so that the long wavelength transverse os- 
cillations of the infinite ID soliton, which cause its snake 
instability, are not present in the spectrum of the 2D 
soliton. 

The occurrence of a coupling between bound and un- 
bound states, which is visible in the avoided crossings in 
the stabilization diagram, is worth stressing. The width 
of the avoided crossings is directly related to the lifetime 
of the bound (resonant) states associated with their de- 
cay into Bogoliubov sound modes. For typical resonant 
states, like the one in Fig. [2j we find a width of the res- 
onance, Acj, of the order of 5 x 10 -3 In [1Q[ we 
noticed that the possible existence of bound states with 
infinite lifetime and dispersion law ~ k 3 ^ 2 could affect 
the thermodynamics of a 2D condensate. The present 
analysis suggests that the bound states of the 2D soliton 
have a finite lifetime, hence making the problem more 
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FIG. 4: Frequency of lowest excited states of the 2D soliton. 
Each point corresponds to an eigenvalue of the matrix equa- 
tion (|13p . whose eigenfunction ip is localized in the soliton 
region and exhibits oscillations along y with wave vector k. 
The solid line is the dispersion law of the stable branch of 
excitations of a ID gray soliton in two dimensions, Eq. J8}. 
All quantities are dimensionless. 

complex. 

In conclusion, a rarefaction pulse is an interesting soli- 
tonic excitation of a 2D condensate, which moves at ve- 
locity close to the Bogoliubov sound and is the low energy 
counterpart of a self-propelled vortex- ant ivortex pair. Its 
momentum and energy were calculated in Ref. [5]. The 
stability of the 2D rarefaction pulse in the KP limit was 
derived in [l5[ from the boundness of the Hamiltonian. 
It was also investigated in [6] , by means of both analytic 
arguments and numerical simulations. More recent nu- 
merical results on the linear stability have been reported 
in However, a detailed calculation of the spectrum of 
the 2D condensate in the presence of the soliton has not 
yet been performed. In this work we have presented the 
results of such a calculation, including an analysis of res- 
onant states based on a stabilization method. We have 
found localized states which closely resemble the stable 
branch of excitations of a ID gray soliton. In the stabi- 
lization diagram, these states appear as resonant states 
coupled to the bath of unbound Bogoliubov phonons. We 
think that these results can be of interest for the inves- 
tigation of ultracold bosonic gases in disk-shaped confin- 
ing potentials, where vortex pairs are known to play a 
crucial role [18[ . The observation of solitons in these sys- 
tems would represent a nice manifestation of nonlinear 
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dynamics in low dimensional superfluids. 
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